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ABSTRACT 

The Galactic globular cluster NGC 1851 has raised much interest since HST photometry re- 
vealed that it hosts a double subgiant branch. Here we report on our homogeneous study into the 
cyanogen (CN) bandstrengths in the RGB population (17 stars) and AGB population (21 stars) using 
AAOmcga/2dF spectra with R ~ 3000. We discover that NGC 1851 hosts a quadrimodal distribution 
of CN bandstrengths in its RGB and AGB populations. This result supports the merger formation 
scenario proposed for this cluster, such that the CN quadrimodality could be explained by the super- 
position of two 'normal' bimodal populations. A small sample overlap with an abundance catalogue 
allowed us to tentatively explore the relationship between our CN populations and a range of elemen- 
tal abundances. We found a striking correlation between CN and [O/Na]. We also found that the 
four CN peaks may be paired - the two CN-weaker populations being associated with low Ba and 
the two CN-stronger populations with high Ba. If true then s-process abundances would be a good 
diagnostic for disentangling the two original clusters in the merger scenario. More observations are 
needed to confirm the quadrimodality, and also the relationship between the subpopulations. We also 
report CN results for NGC 288 as a comparison. Our relatively large samples of AGB stars show that 
both clusters have a bias towards CN-weak AGB populations. 

Subject headings: Stars: AGB and post-AGB — globular clusters: general — globular clusters: indi- 
vidual(NGC 1851, NGC 288) 



1. INTRODUCTION 

Galactic globular clusters (GCs) are no longer thought 
to be perfectly homogeneous, simple stellar populations. 
Although almost all are chemically homogeneous with 
respect to Fe and heavier elements (omega Cen, M22, 
Terzan 5 and NGC 1851 being exceptions), it has long 
been known that GCs show large star-to-star abundance 
variations for l ight elements (e.g., C, N, O, Na. Mg, A l; 
see reviews bv lKraftJll99l IGratton eFail 12004. 2012al). 
These inhomogeneities are considered anomalous because 
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they are s een in very few halo field stars of simila r 
metallicity (jGratton et al.ll2000tlMa,rtell fe Grebell2010l) . 
Studies of all phases of evolution, including the red gi- 
ant branch (RGB), main sequenc e (MS) and subgian t 
branch (SGB; e.g. in NGC 6752. IGratton et alJ[200l . 
have shown the same anomalies. This suggests that many 
of the abundance variations arose in the early phases of 
cluster evolution. 

Recently it was discovered that the globular clus- 
ter NGC 1851 has a double SGB, whereby two evolu- 
tionary sequences are c learly visible (HST photometry, 
iMilone et ai]|2008ll2009t) . The RGB has also been shown 
to split into two when using particular filters ([Han et al.l 
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2009). Spectroscopic observations o f this cluster show 
bimodality in s-process abundances (lYong fc Grundahll 
[20081 iVilianova et~aT1 [20101 iGratton et al.l l2012bO and 
a small spread in [F e/H] (rms scatter ~ 0.05 dex, 
iCarretta et~aD l2011h . IGratton et"aTll2012bl report that 
the two SGB populations ha ve slightly different h eavy 
element contents on average. ICarretta et "all (|201lD find 
they can split their RGB sample into a metal-rich and 
metal-poor population based on the Fe-Ba plane, and 
that each population has its own O-Na anticorrelation. 
This ties in well with the horizo ntal branch (HB) obser- 
vations of lGratton et al.l (|2012d ) who also report two sep- 
arate O-Na anticorrelations. As an explanation for these 
abundance anomalies, as well as the bimodal HB, it has 
been suggested that NGC 1851 may be a product of a 
merg er between two GCs (Ivan den Berghlll996t iCatelanl 
[19971 ICarretta et all 120101 ). iBekki fc Yond ()2012t ) re- 
cently showed that a merger scenario for NGC 1851 is 
dynamically plausible. For reviews on t he phen o meno n 
of mult i ple po pula tions in GCs s e e e.g. iPiottol (|2009f ). 
iMartelll (f201lh and lGratton"et~aTI (j2012al) . 

One of the first inhomogeneities discovered in glob- 
ular clusters was that of the molecule cyanogen (CN, 
often used as a proxy for nitrogen). A pi cture of 'CN- 
bimodality' emerg e d in th e 1970s and 80s (lHesserlll978l; 
INorris fc Freernanl fl979l: ICottrell fc Da Costal Il98lfl 
whereby stars in one population show weak absorption 
by CN ('CN-weak' stars) and stars in the the other show 
strong absorption by CN ('CN-strong' stars). This has 
been observed in most, if not all, clusters. With the 
recent interest in NGC 1851 there have been a c ouple 
of studies of C N on the MS dPancino et al.l I20T0I) and 
the two SGBs (|Lardo et all 12012ft . There does how- 
ever appear to be a dearth of studies of CN in giants 
in NGC 1851 - here we report on observations focusing 
on cyanogen band strengths in the RGB and AGB stars 
of NGC 1851. 

2. STELLAR SAMPLE, OBSERVATIONS AND DATA 
REDUCTION 

Our stellar sam ple was taken from the BV photom- 
etry catalogue of [Walkfi] Cl992|). This catalogue was 
chosen because the photometry is precise enough to dis- 
tinguish between the RGB and AGB populations and 
because it provides accurate astrometry, an important 
feature for multi-fibre spectroscopy. The chosen sample 
of stars was cross-ma tched with the 2MASS catalogue 
(jSkrutskie et al.ll2006f ). thus all of our sample stars are 
actually 2MASS objects, with positions accurate to ~ 0.2 
arcsec. Since the two giant branches merge in the CMD 
at higher luminosities we limited our RGB and AGB sam- 
ples to V > 14 .2. We show our program stars against the 
Walkej] (| 19921 ) CMD in Figure Q] and provide a list in Ta- 
ble [TJ 

Our observational data were taken over the sec- 
ond halves of the nights of the 5th, 7th, 8th and 
9th of September 2009 at the A AT using the multi- 
fibre spectrograph, AAOmega/2dF dLewis et al.1 l2002t 
iSaunders etall 12001 iSharp et all |2006|). A total of 9 
hours of exposures were taken, using 3 field plate setups. 
The 1700B grating was used on the blue arm of the spec- 
trograph, which gave a spectral coverage of 3755-4437A 
and includes the violet CN bands around 3850 — 3880A. 
Spectral resolution in this region was R ~ 3000. 
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Figure 1. The s tellar sample. Small dots are all the stars from 
the Walker (1992) CMD. Filled triangles (red) are our sample of 
17 RGB stars, filled squares (blue) are our 21 AGB stars. 
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Figure 2. Measured S(3839) CN index versus magnitude for the 
NGC 1851 stars. Filled triangles (red) are RGB stars, open squares 
(grey) are AGB stars. The dashed line shows a least-squares fit to 
the 5 RGB stars at the lower envelope of the distribution. The solid 
line is the same but offset so that the zero level of <5S(3839) (<5CN) 
is coincident with the star with the lowest <5CN. The definition of 
<5CN is shown by the arrow. A characteristic error bar for S(3839) 
is shown at top right (see text for details). 

Data reduction was carried out using the 2dF pipeline 
software 2dFdr (version 3.211, April 2009) provided by 
the AAO. Tram- map fits to the multiple spectra from 
each plate were checked by eye, as were the arc reductions 
and final reduced science spectra. Our final sample of 
spectra contained 17 RGB and 21 AGB stars. 

To quantify the CN band strengt hs in each sp e ctrum 
we used the S(3839) CN index of INorris et all (|l98l 
which compares a section of the CN bands with a neigh- 
boring pseudo-continuum (Eqn. [T]). IRAF was used to 
measure the integrated fluxes of Equation [T] for all the 
program stars. 

p3883 j ,a 

5(3839) = -2.5 log { 3 3 ^ A (1) 

J3883 ^ dX 
3. RESULTS AND DISCUSSION 

In Tabled] we list the S(3839) measurements for each 
star, while in Figure[2]we plot them versus V magnitude. 
CN absorption is known to have a temperature depen- 
dence so we have 'de-trended' the dat a in the same man- 
ner as m previou s CN studies (e.g. INorris et al.l ll98lt 
llvans et all 119991: iMartell et all 120081) by fitting a line 
to the lower envelope of the observations. The value 
(5S(3839) (hereafter 5CN) is then the distance from this 
line to each data point. The resultant SCN values are 
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shown in the lower panel of Figure [3] and listed in Table 
[TJ Errors in wavelength calibration or doppler offsets due 
to velocity dispersion were checked and found to be of or- 
der ~ 10 -3 in SCN. The much larger characteristic error 
bar for S(3839) given in the figures (±0.02) reflects the 
typical differences between measurements of S(3839) in 
two observations of the same star. These pairs of obser- 
vations were taken on different nights and with different 
field plates (other clusters in our broader observational 
campaign were used for this, but the data were taken 
during the same timeframe as for NGC 1851). We found 
that this was by far the largest source of error. This is 
probably to be expected since it reflects the combination 
of many sources of error, including the uncertainties of 
fibre placement, fibre throughput, slight pointing errors, 
seeing variability, as well as errors in the data reduction 
(for example). 

In the upper panel of Figure|3]we show a kernel density 
estimate (KDE) histogram of <5CN based on a Gaussian 
kernel with a bandwidth of 0.035. This bandwidth was 
found to be optimal taking into account the error bars 
and small number statistics. It was tested with many 
other similar datasets (Campbell et al., in prep.). It 
can be seen in Figure El that the distribution of SCN 
in the RGB stars is quadrimodal, having four peaks. 
Changing the KDE bandwidth within reasonable limits 
(~ 0.02 — > 0.05) does not alter this result. This was an 
unexpected result because most clusters in our greater 
sample and in the literature show bimodal distributions. 
As a comparison, in the right-hand panels of Figure [3] we 
show the same plot but for NGC 288. The data for this 
cluster were taken in the same observing run and with 
the same instrument as our NGC 1851 data. We chose 
this cluster for comparison because it has a similar metal- 
licity to NGC 1851 ([Fe/H] 1.3 versus -1.2, respec- 
tively) but a different horizontal branch (HB) morphol- 
ogy - NGC 1851 has a bimodal HB morphology while 
NGC 288 has a blue HB only. As can be seen in this 
Figure we find the standard bimodality in SCN in the 
RGB population of NGC 288. 

There is the possibility that the observed quadrimodal- 
ity in our sample of NGC 1851 RGB stars is a chance 
occurrence due to a small sample size. A small random 
sample drawn from a bimodal distribution may give this 
result. To estimate the likelihood of this happening we 
conducted the following test. We established an arbi- 
trary SCN distribution consisting of two Gaussians with 
centres at 0.05 and 0.45 and FWHMs of 0.10. Samples 
of 17 stars were then randomly drawn from this distri- 
bution and the pseudo data was then smoothed with a 
Gaussian kernel having a bandwidth of 0.035. We re- 
peated this 1000 times. The results from this test in- 
dicate that a small fraction (~ 2.5%) of the randomly 
generated samples do present quadrimodal distributions 
(although only 1 out of 1000 was as clearly defined as in 
the real data). Therefore we cannot eliminate the small 
probability that the observed SCN distribution in the 
RGB population has been drawn from a bimodal popula- 
tion through chance. However the case for quadrimodal- 
ity is strengthened by the AGB results in Figured]- the 
AGB population also appears to be quadrimodal in SCN. 
Here the the stars are distributed differently though, with 
the majority of the AGB stars being in the two CN- 
weaker subpopulations, as is usually the case in GCs (e.g. 



NGC 288 in Fig. 0} INorris et al.l[l98lt ICampbell etlll 
1201 If ). The probability for attaining a quadrimodal dis- 
tribution in the RGB population and simultaneously in 
the AGB population by chance then becomes extremely 
small, since they are essentially independent populations 
with different internal distributions of <5CN. 

It has been suggested that NGC 1851 may be a merger 
product between two GCs, initially as an explanation for 
produ cing its bimodal HB (jvan den B ergh 199l lCatelaf] 
fl997h . In this scenario it would also be expected that 
each merging population would have two 'normal' sub- 
populations, each with its own C-N and O-Na (and pos- 
sibly Mg-Al) anticorrelations, and that the superimposi- 
tion of these populations woul d present dual a nticor re- 
lations. In the RGB study by ICarretta etaLl (|2011f ) it 
is indeed found that NGC 1851 has two Na-0 anticor- 
relations, one in their metal-rich population and one in 
their metal-po or population. This t ies in well with the 
findings of the iGratton et al.l (|2012d ) study where it was 
found that there are also two independent Na-0 anti- 
correlations on the horizontal branch - one in the RHB 
population and one in the BHB populatio n. With re- 
gards to C and N the picture is less clear. lLardo et al.l 
(2012) studied the SGB populations and found a spread 
in C and N between stars. These elements were also 
found to be anticorrelated, however no bimodal signa- 
ture was detected. The resolution of the spectra in that 
study was ho wever quite low , with R ~ 1000. In the 
SGB study of IGratton et all (|2012bT) it was found that 
the two subgiant branches have different average C abun- 
dances. Interestingly they also found that there are dif- 
ferent proportions of C-normal and C-poor stars in each 
SGB, which may indicate that each SGB hosts multiple 
subpopulations, again suggestive of a merger scenario. 
In a s tudy of the two RGB populations Vill anova et al.l 
(2010) found a spread in CNO elements but f ound no dif- 
ferenc e in C+N+O between the populations. lYong et ail 
(2009) did however find a significant variation (a factor 
of 4) in C+N+O in four RGB stars. Thus there is still 
uncertainty as to whether C+N+O is constant between 
populations or not. This is an important diagnostic since 
it is a very useful discriminant between possible polluters 
of the primordial material from which N-rich populations 
form. Cluste r age determinations a r e also very sensitiv e 
to C+N+O (|Rood fc CrockeHfl985l: iCassisi et al.l [20081) . 
More information/observations of CNO are needed to 
clarify the situation. 

In the case of CN bandstrengths the merger scenario 
leads to a natural expectation that the two bimodal 
CN populations in the original clusters would also su- 
perimpose, giving a quadrimodal distribution. Thus 
our discovery of a quadrimodal distribution of CN band 
strengths in the RGB and AGB populations of NGC 1851 
adds further weight to the merger formation scenario for 
this cluster. The CN quadrimodality also suggests that 
there may be four populations with different N abun- 
dances in this cluster. Although CN is generally accepted 
as a proxy for N we note that the band strengths may 
also be affected by the abundances of C and O. Thus, 
again, a complete set of (absolute) abundance observa- 
tions including C, N, and O are needed. 

With a view to gaining more information about the 
four SCN popul ations we performed a cross-match of our 
sample with the ICarretta et all (|201lD catalogue of RGB 
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Figure 3. Left-hand Bottom Panel: <5S(3839) CN index versus V magnitude for NGC 1851 (see Fig. [2] for the definition of <5S(3839)). 
Left-hand Upper Panel: Normalised kernel density estimate (KDE) histogram (Gaussian kernel, bandwidth = 0.035) of 5S(3839). A clear 
CN quadrimodality can be seen on the RGB and the AGB. A characteristic error bar for S(3839) is shown at top right of the lower panel. 
Right-hand Panels: Same as left-hand panels but for NGC 288, for comparison. The bandwidth for the Gaussian kernel is the same as 
that used for th e NGC 1851 data. In this case photometry (uvby) for the stellar sample selection was provided by F. Grundahl (private 
communication; Grundahl ct al. 1999) 



Figur e 4. NGC 1851 RGB stars in common with the 
Carrctta ct al. (2011) catalogue. Various abundance ratios plus 
radial velocity from the catalogue are shown versus our <5S(3839) 
values. The vertical dotted lines denote the peaks of the 4 5CN 
populations (see Fig. Q5). The horizontal lines in [Fe/H] and radial 
velocity (RV ) show the cluster m eans of —1.16 and +320.26 km/s 
respectively ([Carrctta ct al. 2011]). 



star abundances. The catalogue contains 124 stars with 
a range of abundance measurements for p-capture, a- 
capture, Fe-peak, and n-capture elements. Our cross- 
match found an overlap of 10 RGB stars and one AGB 



star (Table [[]). Not all of the ICarretta et~aH (|2011[ ) 
objects have the full set of abundance results but for 
some abundances we have 7-10 stars to compare with. 
In Figure [4] we sho w various abundance ratios from 
ICarretta etall (|2011[ ) against our RGB <5CN values. In 
the second panel of Figure H] it can be seen that there ap- 
pears to be no correlation between our four <5CN groups 
and [Fe/H]. This is con trary to what we mig ht expect 
from the SGB results of iGratton et all (|2012bl) where it 
was found that the two SGB populations differ in average 
[Fc/H]. We note however that the cross- matched sample 
is very small, especially within each of the four subpopu- 
lations, which contain only 1-4 stars each. Sodium on the 
other hand shows a definite correlation with SCN. This is 
typical of GC abundance anomalies, where N is higher in 
the stars with high Na. Barium also appears to correlate 
with (5CN. In contrast oxygen shows little variation with 
SCN. By using the ratio [O/Na] the 'noise' of the Fe scat- 
ter can be removed, and the relationship between O and 
Na is amplified. We show [O/Na] in the top panel of Fig- 
ure |3J Here there is a striking anticorrelation, such that 
the CN-weak populations show much higher O /Na ratios 
than the CN-strong populations. The lack of correlation 
between Fe (and Ca, Si, Ti - not shown) and the light el- 
ements + s-process elements suggests that the Fe-group 
(+ a capture) nucleosynthetic source(s) are separate to 
the light element + s-process source(s). AGB stars are 
the most likely primordial source for the O, Na, CN and 
s-process enhancements. If the four populations are in- 
deed real, then this suggests that each population was 
polluted by AGB star ejecta to differing degrees. We 
note that the star with intermediate [Ba/Fe] ~ 0.5 is 
the only star in the sam ple with no K magnitude from 
the ICarretta et aT] (|2011l) database and thus may be an 
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Table 1 

List of target stars for NGC 1851. IDs (column 2), V magnitudes and B-V values are from|Walkcr (1992). Column 4 IDs are from 
Carrctta ct al. (2011) and show the overlap between studies. In column 7 are our raw CN band strength measurements, S(3839), and in 

column 8 our de-trended CN index values <5S(3839). 



TYPE 


ID (Walk.) 


ID(2MASS) 


ID(Carr.) 


v 


B-V 


Sl'3839') 


r5Sf38391 


RGB 


24 


(151 34832-40031 51 




15.40 


1.00 


0.31 


0.19 


RGB 


28 


051 34897-4001 1 99 


41113 


15.71 


0.94 


0.34 


0.26 


RGB 


71 


051 3541 4-4003038 


32112 


15.45 


0.98 


0.73 


0.62 


RGB 


79 


051 35462-4005094 




15.29 


1.00 


0.15 


0.01 


RGB 


120 


051 35671 4001 01 6 




15.00 


1.06 


0.57 


0.40 


RGB 


151 


051 35898-395Q586 




16.19 


0.90 


0.43 


0.41 


RGB 


160 


051 35866-40001 78 




16.14 


0.89 


0.42 


0.40 


RGB 


161 


(151 35867 40001 90 


44803 


14.91 


1.07 


0.39 


0.20 


RGB 


162 


05135869 3Q5Q949 


46228 


16 19 


0.91 


0.11 


09 


RGB 


208 


051 35Q77 4001 374 




15.66 


0.93 


0.42 


33 


RGB 


368 


0514095Q 4000990 


44414 


15.62 


0.95 


0.33 


0.24 


RGB 


441 


05140365 4001 5Q6 




16.18 


0.89 


0.28 


26 


RGB 


1028 


051 41 059-3Q580Q5 


47385 


15.41 


1.00 


0.12 


0.00 


RGB 


1256 


051 41 724-4002080 


37070 


14.38 


1.21 


0.34 


0.08 


RGB 


1284 


05141Q56 4004055 


26801 


14.39 


1.18 


0.63 


37 


RGB 


1286 


051 41 947-4000076 


44939 


14.95 


1.09 


0.24 


0.06 


RGB 


1323 


05149981 4001551 


38215 


16.19 


0.89 


0.58 


0.56 


AGB 


182 


05135Q18 40094Q6 




14.93 


0.92 


0.35 


17 


AGB 


222 


051 4001 Q 40099Q1 




14.69 


0.98 


0.39 


0.17 


AGB 


245 


05140068 400393Q 
ujiiuuuo "iuwo^oy 


30315 


14.53 


1.08 


0.24 


0.00 


AGB 


430 


05140355 40094QQ 




15.03 


0.95 


0.16 


-0 01 


AGB 


506 


051 40446 40031 1 3 




14.76 


0.88 


0.29 


0.08 


AGB 


572 


05140508 4009978 




15.56 


0.82 


0.31 


0.21 


AGB 


633 


05140584 4009196 




14.25 


1.16 


0.26 


-0 01 


AGB 


680 


05140659-4002026 




14.78 


1.04 


0.44 


0.24 


AGB 


697 


05140701-4003449 




14.73 


1.02 


0.17 


-0.04 


AGB 


711 


05140758-4003164 




15.50 


0.83 


0.22 


0.11 


AGB 


848 


05140883-4002380 




14.23 


1.15 


0.23 


-0.05 


AGB 


849 


05140900-4004539 




14.83 


0.97 


0.41 


0.21 


AGB 


887 


05140916-4002296 




14.99 


0.87 


0.27 


0.09 


AGB 


988 


05141034-4004235 




14.76 


0.98 


0.70 


0.49 


AGB 


989 


05141026-4003150 




15.42 


0.80 


0.17 


0.05 


AGB 


1004 


05141057-4003308 




14.93 


0.94 


0.33 


0.15 


AGB 


1014 


05141035-3958148 




14.72 


1.04 


0.57 


0.36 


AGB 


1037 


05141084-4001475 




15.19 


0.88 


0.52 


0.37 


AGB 


1172 


05141351-4003408 




14.87 


0.96 


0.66 


0.47 


AGB 


1214 


05141501-4003040 




14.85 


0.97 


0.39 


0.20 


AGB 


1246 


05141641-4002214 




15.42 


0.80 


0.29 


0.17 



unreliable data point. If so, then Ba would present a bi- 
modal distribution (panel 3 of Fig. [4j. Importantly each 
mode of Ba abundance would be associated with two 
<5CN peaks: the two CN-weaker populations would be 
Ba-poor compared to the two JCN-richcr populations. 
This could be a useful diagnostic for disentangling the 
multiple populations in the merger scenario. Clear bi- 
modalities in s -process abundances have been reported 
for t he RGB (lYong fc Grundahll [200l iVillanova et al.1 
20101) and the SGBs (iGratton et al.ll2012b| ). Interestingly 
Carretta et al.1 (|2011l ) do not report a bimodality in Ba 
but do show that s-process elements are (anti) correlated 
with p-capture elements. In the bottom panel of Fig- 
ure |4] we show the radial velocities for the cross- matched 
sample. It can be seen that only the most extreme <5CN 
population stands out, having an average radial veloc- 
ity ~ 5 km/s lower than the other three populations, 
which might suggest this group is kinematically distinct. 
Again we stress that this is a very small dataset so the 
discussion above is only speculative. We note that our 
group is in the process of collecting medium resolution, 
broad wavelength coverage s pectra using 2dF/AA omega 
to complement the excellent iCarretta et al.l (|201lD RGB 
dataset with C and N abundances. When complete the 
combined dataset will allow a 'holistic' analysis (includ- 
ing absolute abundances of C, N, O and therefore the 



sum C+N+O) of the abundance and population trends 
for NGC 1851 red giants. 

4. SUMMARY AND CONCLUSIONS 

We have recorded a homogeneous set of spectra for 
17 RGB and 21 AGB stars in the globular cluster NGC 
1851. We find that the CN band strengths divide into 
four groups in both the RGB and AGB populations. This 
lends support to the theory that NGC 1851 formed from 
a merger of two clusters since one of the expected signa- 
tures of this would be two superimposed bimodal distri- 
butions in CN. 

We cross-mat ched our sample with t hat of the high res- 
olution study of ICarretta et al.l (|2011l ) and found a small 
number of stars in common. This gave us the opportu- 
nity to compare elements that typically (anti)correlate 
with N in globular clusters, such as O, Na and Ba. We 
found that Na did indeed correlate with <5CN. An anti- 
correlation [O/Fe] was less clear but when considering 
[O/Na] it was found that there was a very strong anti- 
correlation. A possible correlation with Ba was observed. 
The (anti)correlations between these elements and 8CN 
(and thus presumably N and hence C) suggest that the 
material from which each of the four populations formed 
was polluted by AGB stars. We also speculated that the 
Ba distribution may be bimodal, as found in previous 
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studies. If so, then the two CN-weaker and two CN- 
richer populations would be paired, and this may reflect 
a distinction between the two GCs in the merger hy- 
pothesis. It must be noted that the comparison sample 
is small, so strong conclusions could not be made. Large- 
sample, high-resolution observations combining absolute 
abundances of C, N, O, Fe, neutron-capture elements and 
radial velocities in the giant branches of NGC 1851 are 
needed to check that there are indeed four chemically 
(or even kinematically) distinct populations and to de- 
termine which sub-populations are related to each other. 
Our group is in the process of collecting observational 
data to this end. 

Finally we note that the AGB samples we presented 
here are the largest AGB samples in any GCs to date. 
Our finding that NGC 1851 and NGC 288 both have 
CN-weak dominated AGB populations adds to a grow- 
ing picture in the literature that the AGB CN dis- 
tributions in__GCs__are different to t he RGB distri- 
butions dCampbell et all [2fM [20ll lLai et all l20Tll; 



iSmolinski et all 120111: iSimpson et all 120121 ) For brief 
reviews on thi s topi c see iSneden etaiT (f2000h and 
ICampbell et al.l ([200 6). The current study forms part 
of a larger study with relatively large datasets of AGB 
stars to confirm this for a range of GCs (Campbell et al., 
in prep.). 

SWC acknowledges support from the Australian Re- 
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